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For gravitational-wave interferometric detectors, there is a tradeoff between the detector bandwidth and peak 
sensitivity when focusing on the shot noise level. This has to do with the frequency-dependent propagation 
phase lag (positive dispersion) of the signal. We consider embedding an active unstable filter—a cavity-assisted 
optomechanical device operating in the instability regime—inside the interferometer to compensate the phase, 
and using feedback control to stabilize the entire system. We show that this scheme in principle can enhance the 
bandwidth without sacrificing the peak sensitivity. However, there is one practical difficulty for implementing 
it due to the thermal fluctuation of the mechanical oscillator in the optomechanical filter, which puts a very 
stringent requirement on the environmental temperature and the mechanical quality factor. 


Introductions .—Advanced gravitational-wave (GW) detec¬ 
tors, including Advanced LIGO ifTl. Advanced VIRGO 0 
andKAGRAll, are dual-recycled Michelson interferometers 
with both power-recycling mirror (PRM) and signal-recycling 
mirror (SRM), as shown schematically in Fig.[TJa). The GW- 
induced differential motion of the two end test masses (ETMs) 
modulates the laser at frequency (Oq and creates sidebands at 
(Oo ± fi with fi being the GW frequency, ranging from lOHz 
to lO^Hz for ground-based detectors. In the presence of SRM, 
these sidebands are coherently reflected back to the interfer¬ 
ometer, which modifies the detector frequency response and 
bandwidth. For instance, advanced LIGO, in its nominal oper¬ 
ation mode, is using SRM to broaden the detector bandwidth. 
However, this is at a price of decreasing the peak sensitiv¬ 
ity when considering the shot noise, as illustrated in Fig.[2b). 
Such a tradeoff between the bandwidth and the peak sensi¬ 
tivity, known as the Mizuno theorem ||4l, can be attributable 
to the frequency-dependent propagation phase of sidebands 
(positive dispersion), as cavity resonant condition is satisfied 
only for a single frequency. There are several approaches pro¬ 
posed in the literature aiming at broadening the detector band¬ 
width without degrading its peak sensitivity. One is the so- 
called white-light-cavity idea lEHni that uses atomic medium 
with negative dispersion to cancel the positive dispersion. The 
other is applying either external squeezing uniiii or internal 
squeezing ifTSl [Tbl to reduce the shot-noise level while keep¬ 
ing a broad detector bandwidth. In particular, external squeez¬ 
ing has already been implemented in large-scale GW detec¬ 
tors iniiisi and also planned for future upgrades nMn. 

We consider a different approach to realize such a broad¬ 
band enhancement by placing an optomechanical device as a 
filter inside the signal-recycling cavity, illustrated schemati¬ 
cally in Fig.[TJa). Such an optomechanical filter operates in 
the unstable regime pumped by an additional laser, and the 
entire scheme is stabilized via feedback control. Within the 
signal loop, the unstable filter acts as a phase compensator 
with negative dispersion. Notice that there is no violation of 
causality, as the controlled system is stable and the measured 
output always lags behind incoming GW signals. In addition, 
the feedback uses the SRM output, as indicated by the flow 
chart in Fig.[TJc), which contains the GW signal together with 
noise, and therefore it does not influence the signal-to-noise 


ratio, as proven in Refs. ll22l l23l . This idea is inspired by 
recent experimental observations of the optomechanical ana¬ 
logue of the electromagnetically induced transparency by S. 
Weis et al. Il24l . Teufel et al. Il25l and Safavi-Naeini et al. 1^ . 
and also a more recent theoretical proposal by Ma et al. Wl- 
These are, instead, using the optomechanical interaction to en¬ 
hance the phase lag within a certain frequency range. 

A brief summary .—Before develling into the details, we 
summarize the main features and explain them qualitatively. 
For the optomechanical filter, the radiation pressure couples 
the optical field intensity and the mechanical displacement. 
Such a nonlinear coupling is analogous to the three-wave mix¬ 
ing in nonlinear optics. By tuning the pump laser frequency 
to be fi)o + (Om with the mechanical oscillator frequency (Om 
much larger than the cavity bandwidth jf, i.e., 

( 1 ) 

the cavity resonance at Oo is in favor of the down-conversion 
process, which amplifies the sidebands around (Oq and me¬ 
chanical motion. It can be viewed as a phase-insensitive para- 
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FIG. 1: (a) The dual-recycled Michelson configuration for advanced 
gravitational-wave detectors, with an unstable optomechanical filter 
(blue) embedded. An internal SRM (iSRM) is introduced to match 
the transmissibility of input test mass (ITM) mirror to remove the 
arm cavity pole; (b) The tradeoff between the detector bandwidth 
and peak sensitivity (black and red curves), and the effect of adding 
unstable filter in the ideal scenario (blue); (c) A simple flow chart. 
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metric amplifier for sidebands with the following input-output 
relation; 


^out(^) 


+ + 7opt) 

^ ii^Ym Tipt) 


Qin(^) ■ 


( 2 ) 


Here 7 „, = (Om/Qm is the mechanical damping rate with Qm 
being the quality factor; 7opt is the negative mechanical damp¬ 
ing rate due to the optomechanical interaction, and is approxi¬ 
mately equal to PcJ^ (Oq / {m( 0 ,„c^), where Pc is the intra-cavity 
laser power, ^ is the cavity finesse, and m is the mass of the 
mechanical oscillator. 

The unstable regime we referred to is when 7opt becomes 
much larger than 7 „, and the mechanical damping rate be¬ 
comes negative. With the feedback control engaged, the open- 
loop input-output relation of the optomechanical filter is 


flout(ii) « ^'^'.^°P‘ flin(n) « - exp — ^ flin(ii), (3) 
i2-i7opt V 7opt/ 


which exhibits negative dispersion rather than positive disper¬ 
sion for a passive filter cavity. To compensate the propagation 
phase lag 0ann = 2nLann/c with Larin being the arm cavity 
length and equal to 4km for LIGO, we therefore require 

7opt — cj Larm , (4) 


and the resulting intra-cavity laser power scales as: 


Pc 


\Tarm/ \ t^/ 27 r / \ m J 



(5) 


In Fig.|^ we show the resulting shot-noise-only sensitivity 
with nominal parameters given above and one additional pa¬ 
rameter Lf — 1cm being the filter cavity length. We can see 
that the unstable filter can compensate the phase at low fre¬ 
quencies but not perfect at high frequencies due to high-order 
frequency dependence of the phase delay, which in principle 
can be improved by cascading several unstable hlters. 

So far, we have only mentioned the phase property of the 
unstable filter. As a general principle proven by Caves for 
phase-insensitivity parametric amplihers ll28l . there will be an 
additional noise term given Eq. (|^; 


iladd(^) — 


2.^/7/! 7)pt 

^ T ^{Ym 7opt) 




( 6 ) 


As shown later, it comes from the coupling of the mechani¬ 
cal oscillator to the environmental thermal bath. The spectral 
density for Lth is approximately given by -I- 1 

with Lenvir being the temperature. In order for the thermal 
noise to be lower than the quantum shot noise, we require 


SkBLenvir /Qm < ^7SRM , (7) 

with 7srm = cLsRM/(4Larm) and Lsrm being the power trans¬ 
missivity of the SRM. Here 7 srm is the original detector band¬ 
width, as ITM and iSRM are impedance matched and the 
bandwidth is solely determined by SRM. Therefore, the lower 
is the detector bandwidth that we start off, the higher require¬ 
ment will be imposed on the temperature and quality factor. 


■> 


c 

u 


c3 


ca 

si 

PLh 



Frequency (Hz) 


FIG. 2: The top figure shows the shot-noise-only strain sensitivity 
of the unstable-filter scheme, compared with the broadband and nar¬ 
rowband tuned cases without the filter. The power and arm cavity 
length are the same as Advanced LIGO. At high frequencies, its sen¬ 
sitivity deviates from the ideal scenario because of imperfect phase 
cancellation of higher-order terms of £1 in (bottom figure). 


As an order of magnitude estimation, we have 


Tenvir 

Qm 


<6x 


7sRM/27r \ 
lOOHz ) ■ 


( 8 ) 


In Fig.[^ we show the thermal noise effect on the sensitivity, 
which is significant. As mentioned in Ref. If27l . one approach 
for mitigation is applying the optical-dilution idea Il29l432l . 
which uses optically-induced rigidity to dilute the mechanical 
dissipation, and allows for enhancement of Q,„ by a factor of 
hundred or even more but further experimental study is neces¬ 
sary. After this summary, we will present more details about 
the scheme and also the issue of feedback control. 

Dynamics of the optomechanical filter .—We start with the 
optomechanical filter shown in Fig.j^a), derive the input- 
output relation for the sidebands, and later combine it with 




FIG. 3: Effect of the thermal fluctuation of the mechanical oscillator 
on the sensitivity. 



























3 


(a) 


I—I ^0 I— 

flout ^ ^ _ 


(b) 


UJ'jyi 


C^o 


CJO + OJrt 


FIG. 4: (a) schematics of the optomechanical filter; (h) frequency of 
interests: cavity resonance (Hq (red), laser frequency COQ + Om (blue) 
and its lower sideband (black), with the upper sideband around Wq + 
2cOm not shown. 
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the main interferometer. It is an optomechanical device which 
has been studied extensively in the literature (see recent re¬ 
views: Refs. 133] El), with Hamiltonian given by — 
Ho + Hint + Hy^. -t- Hy^. The free part Hq is 

„ r P ( 9 9 

Ho = hcooa'a + - - \--mco^x , (9) 

2m 2 


The linearized interaction Hamiltonian Hint is 


Hint = -hgo 


n^i((Bo+(l),„)r _j_ 


( 10 ) 


with go = (Ooa/Lf and a = (2PcLf/(hcOoc)y^^; Hy^ describes 
how the cavity mode d interacts with ingoing (outgoing) field 
flin (flout); Hy„ describes the coupling between the mechanical 
oscillator with the environmental thermal bath. 

The parameter regime we are interested in, as mentioned 
earlier in Eq. ([T]), is the so-called resolved-sideband regime 
illustrated in Fig.|l^b). This allows us to ignore the upper me¬ 
chanical sideband around (Oq -F 2®,,, and use the rotating-wave 
approximation (RWA) in the interaction picture, obtaining 

HM^ = -hgidb+dny, (11) 

where we have introduced the annihilation operator b for the 
mechanical oscillator through x{t) = Xq[be^‘“’'”‘ +b^e‘“”'‘) 
and g = goXq with Xq being the ground-state uncertainty. 

The resulting Heisenberg equations of motion read: 

d(t} + yfd(t) = igb'^(t) +y^din(t), (12) 

b(t) + Y„,b(t) = igdyt) + ^/^bth(t). (13) 

Solving them in the frequency domain, we obtain the input- 
output relation for sidebands using dout = ~flin + \/2Yfd: 

flout(n) « ^^y"^^°P‘j fl^in(n)-F«add(n) , (14) 

where we have used Yf ^ ^ and defined 7opt = g^/Yf ~ 

Pc^(Oo/{m(OmC^)- The first term gives Eq. (|^ and in Eig.|^ 
we compare it with the exact phase of the optomechanical fil¬ 
ter without RWA and Yf ^ The second term finii is the 
thermal noise term mentioned earlier. 

Eor the mechanical sideband, the resulting susceptibility is: 

~ (ifi-f Ytn Tipt) • (15) 


FIG. 5: The difference between obtained without using RWA 
and —2arctan(f2/y(jpt). The parameter regime with co„, ^ Yf^ 7opt 
is thus preferred for matching the required frequency dependence. 


The mechanical motion is unstable in the parameter regime of 
interest with 7opt much larger than the bare mechanical damp¬ 
ing rate Ym, and a feedback control is thus needed. We show 
the control scheme after considering the filter together with 
the main interferometer, as given below. 

Dynamics of the entire system .—The total Hamiltonian 
reads Htot = ^ifo + ^filter + ^ifo-fiiter- We can model the main 
interferometer also as an optomechanical device lESlI : 

^ ^ P^ 

Hifg = hcood^dHy. H- hGoid-\-d'')X ~\-XpQyy. (16) 

Here d is the differential optical mode —a single-mode ap¬ 
proximation which is valid when considering sideband fre¬ 
quency much lower than one free spectral range c/(2Lann); 

describes the interaction between d and the ingoing (out¬ 
going) field din{dout) at the dark port; Gq = Cdotf/Tarm 1^ 
the coupling strength and d = (2f’annTarm/(^®oc))'/^ with 
farm being the arm cavity power; X is the differential mo¬ 
tion of the ETMs and is driven by the GW as a tidal force: 
Fgw = MLnmih{t) with h being the GW strain. The coupling 
between d and a can be quantified by (exchanging photons): 

Afo-fiiter = h(oyd^d-Yddy , (17) 

where (Os = cYf /Tarm is the coupling rate, and equal to the 

optomechanical coupling rate g when Eq. 0 is satisfied. 

To focus on the shot noise, we will first ignore the radia¬ 
tion pressure effect on ETMs by assuming M ^ (finite test 
mass is considered in the rigorous treatment). The resulting 
equations of motion under RWA are: 

a = —icOsd -F igb^ , ( 18 ) 

b'^ = -Ymb'^-igd-\-s/2^bly, (19) 

d = —7SRMt( — i(Osd -F \/27sRMt(m + iGoLnrmh . (20) 

The input-output relation at the SRM is liout = <7in — v/27srm d. 

The system stability can be examined from the eigenvalues 
of the dynamical matrix (read off from equations of motion); 


A = 


0 ig -iCOs 
-ig -Ym 0 
—/®v 0 — 7 srm 


( 21 ) 

























4 


Having eigenvalues with positive real part implies instability, 
which is the case given the relevant parameter regime. 

To find the stabilizing controller, we follow the state-space 
approach. We first need examining observability and control¬ 
lability of the system. It turns out that if we detect the output 
phase quadrature, which contains the GW signal, using homo¬ 
dyne detection, only one of the two mechanical quadratures 
will be observed, and this will lead to an uncontrollable sys¬ 
tem. An apparent solution is measuring both the amplitude 
and phase quadratures with heterodyne detection. However, 
this is at a price of increasing the shot noise by v/2 in am¬ 
plitude Ofil . Instead, we can pick off a small portion of the 
output signal and use the heterodyne detection only for con¬ 
trol purpose. The rest is still measured using the homodyne 
detection for exacting the GW signal, which will not degrade 
the sensitivity by an noticeable amount. 

The resulting readout vector with heterodyne detection in 
the sideband picture is D = (0, 0, 1)—the measured doui is 
linear to d, and the control input vector is B = (0, 1,0)^ (su¬ 
perscript ^ for transpose)—the feedback force is coupled to 
the mechanical displacement that is linear to fi'. The sys¬ 
tem becomes both observable and controllable, as we have 
rank([D;DA;DA2]) = rank([B,AB,A^B]) = 3. A stabiliz¬ 
ing controller can then be constructed (see section 9.4 of 
Ref Ell), which has the following transfer function (from d 
to b^) in the frequency domain: 

C(f2) = -K(-/m-A-fBK + LD)^'L, (22) 

where K = {Ki,K 2 ,K^) and L = {L\^L 2 tLt,Y are chosen 
such that eigenvalues of A — LD and A — BK all have neg¬ 
ative real part. Given the nominal parameter specification 
in Eq. (|^, the system will be stabilized by setting K = 3 x 
10^e^'(—!, 1, —1) andL = 5 x 10^e^'(!, 1.2,1), where £ is 
the fraction of the output (in amplitude) measured using het¬ 
erodyne detection. With this set of K and L found, arbitrary 
stabilizing controller can be generated via Youla-Kucera pa¬ 
rameterization 1381, which is used for control optimization. 

Rigorous treatment .—We have used several approxima¬ 
tions in order to gain intuitive understanding. These approx¬ 
imations are reasonable in the parameter regime that we are 
focusing on. We briefly outline how the more rigorous treat¬ 
ment is applied in the actual analysis for producing the results 
presented in this paper. Specifically, the optomechanical inter¬ 
action in the main interferometer mixes the upper and lower 


sidebands around coo, which was ignored by assuming infinite 
test mass and focusing on the shot noise only. The interac¬ 
tion in the optomechanical filter mixes in sidebands around 
COq + 2(Om that were ignored in RWA. A rigorous treatment 
therefore involves propagating four sidebands at frequency 
0)0 ± and Oo + 2(0m ± Cl. The result is shown in Fig.|^ as 
one example. Notice that the low-frequency radiation pressure 
noise can be reduced by either using frequency-dependent 
readout ll39l or increasing the test mass size. 

Conclusions .—We have considered using an unstable op¬ 
tomechanical filter to enhance the bandwidth of advanced 
gravitational-wave detectors. The entire system is stabilized 
via a feedback control. Within the loop, phase lag of the side¬ 
bands due to free propagation inside the arm cavity is com¬ 
pensated by the unstable filter. The issue for implementing 
this scheme is the stringent requirement on thermal noise of 
the mechanical oscillator in the filter, and could be mitigated 
by using the optical-dilution idea but requires further experi¬ 
mental verification. 
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FIG. 6: The figure shows the total quantum noise (including radiation 
pressure noise) of the unstable filter scheme at different temperature. 
The color code for different curves is the same as Fig.[^ The test 
mass is assumed to be 40kg, the same as Advanced LIGO. 
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